The Prograde Orbit of Exoplanet TrES-2b 



< 



X 



Joshua N. Winn^, John Asher Johnson^'^, Norio Narita^, Yasushi Suto^, 
Edwin L. Turner^, Debra A. Fischer/, R. Paul Butler^, Steven S. Vogt^ 



5 



Francis T. 0'Donovan^°, B. Scott Gaudi^^ 



00 

O ■ ABSTRACT 
u , 

We monitored the Doppler shift of the GOV star TrES-2 throughout a transit 
of its giant planet. The anomalous Doppler shift due to stellar rotation (the 
Rossiter-McLaughlin effect) is discernible in the data, with a signal-to-noise ratio 
of 2.9, even though the star is a slow rotator. By modeling this effect we find that 
^ . the planet's trajectory across the face of the star is tilted by —9 ± 12 deg relative 

to the projected stellar equator. With 98% confidence, the orbit is prograde. 

O 
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1. Introduction 

A small fraction of Sun-like stars have giant planets with orbital periods smaller than 
about 10 days (Marcy et al. 2005, Udry & Santos 2007). The existence of these planets 
was a surprise, because it was expected that giant planets would only be found beyond the 
"snow line," with orbital distances greater than a few astronomical units. Other surprises 
have come from detailed studies of individual objects. Some are found on highly eccentric 
orbits (Johnson et al. 2006, Bakos et al. 2007, Maness et al. 2007, Johns-KruU et al. 2008). 
Some have mean densities that are quite small (Knutson et al. 2007, Mandushev et al. 2007) 
or large (Sato et al. 2005, Torres et al. 2007) in comparison with Jupiter. 

However, in at least one sense, the close-in giant planets have fulfilled prior expectations: 
they orbit their host stars in the prograde direction, relative to the sense of the stellar 
rotation. This is true, at least, of the 6 systems for which measurements of spin-orbit 
alignment have been reported (Queloz et al. 2000; Wolf et al. 2007; Narita et al. 2007a,b; 
Loeillet et al. 2007; Winn et al. 2005, 2006, 2007a). In all of these cases but one, the sky 
projections of the orbital axis and the stellar rotation axis are observed to be fairly well- 
aligned, with measurement precisions ranging from about 1.5 to 30 deg. The exception is 
HD 17156, for which the angle between those axes was found to be 62 ± 25 deg (Narita 
et al. 2007b). In all of these cases, the measurement technique rehes upon the Rossiter- 
McLaughlin (RM) effect, the anomalous Doppler shift that occurs during transits due to 
stellar rotation (see, e.g., Queloz et al. 2000, Ohta et al. 2005, Gimenez 2006, Gaudi & Winn 
2007, Winn 2007). 

A close alignment between the orbital and rotational axes seems natural because this 
pattern prevails in the Solar system, and because the angular momenta of the parent star 
and the planetary orbits presumably derive from the same protostellar disk. However, 
some theories of planetary migration — proposed to explain how giant planets attain short- 
period orbits — predict occasionally large misalignments (Chatterjee et al. 2007, Fabrycky & 
Tremaine 2007, Wu et al. 2007, Nagasawa et al. 2008). These theories, as well as the gen- 
eral history of surprises in this field, provide motivation to continue measuring exoplanetary 
spin-orbit alignment. 

In this paper, we present a measurement of the RM effect for the transiting exoplanetary 
system TrES-2. This system was discovered by O'Donovan et al. (2006). It consists of a 
planet with a mass of 1.2 Mjup and radius 1.2 i?jup orbiting a GOV star with a period of 
2.5 d (O'Donovan et al. 2006, Holman et al. 2007, Sozzetti et al. 2007). It did not stand out 
as a promising RM target because the star is relatively faint {V — 11.4) and is a slow rotator 
{v sin ii, — 2.0 ± 1.5 km s~^; O'Donovan et al. 2006). On the other hand, the transit occurs at 
a high impact parameter across the stellar disk {b — 0.8540 ± 0.0062; Holman et al. 2007), a 
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favorable circumstance for this type of measurement (Gaudi & Winn 2007). Furthermore, in 
our continuing effort to measure the spin-orbit angles for a statistically meaningful number 
of systems, we do not want to ignore stars with small sky-projected rotation rates. This is 
because a small value of vsini^ might be caused by a small value of sini^, i.e., there might 
be a large spin-orbit misalignment. For these reasons, we pursued TrES-2. We describe the 
new data in § 2, the model that we used to interpret the data in § 3, and the results in § 4. 

2. Observations and Data Reduction 

We observed a transit of TrES-2 on UT 2007 April 26 with the Keck I 10m telescope 
and the High Resolution Echelle Spectrometer (HIRES; Vogt et al. 1994). We set up the 
instrument in the same manner that has been used consistently for the California-Carnegie 
planet search (Butler et al. 1996, 2006). In particular we employed the red cross- disperser 
and used the I2 absorption cell to cahbrate the instrumental response and the wavelength 
scale. The slit width was 0'.'85 and the typical exposure time was 3-4 min, giving a resolution 
of about 70,000 and a signal-to- noise ratio (SNR) of approximately 200 pixel"^. We observed 
the star for 4 hr bracketing the predicted transit midpoint and obtained a total of 56 spectra, 
of which 30 were taken during the transit. 

We also obtained two iodine-free spectra, with a higher SNR and higher resolution. We 
used the sum of these spectra as a template for the Doppler analysis, which was performed 
with the algorithm of Butler et al. (1996). We estimated the measurement error in the 
Doppler shift derived from a given spectrum based on the scatter among the solutions for 
individual 2 A sections of the spectrum. The typical error was 6 m s~^. The data are given in 
Table 1 and plotted in Figs. 1 and 2. Also shown in those figures are data obtained previously 
by O'Donovan et al. (2006), consisting of 11 velocities measured with Keck/HIRES using a 
different setup0, as well as the photometric data of Holman et al. (2007). 

3. The Model 

To determine the projected spin-orbit angle and its uncertainty, we simultaneously fitted 
a parametric model to the radial- velocity data as well as the photometric data of Holman et 
al. (2007). We included the photometric data as a convenient way to account for the uncer- 



^Table 3 of O'Donovan et al. (2006) gives incorrect values for the heliocentric Julian dates of the velocity 
measurements. The corrected dates were provided to us by D. Charbonneau (private communication, 2007). 
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tainties in the photometric parameters and their covariances with the spin-orbit parameters, 
although in practice the photometric uncertainties were irrelevant for this system. 

The model is based on a circular orbit of a star and planet. The photometric transit 
model was identical to the model used by Holman et al. (2007). To calculate the anomalous 
Doppler shift as a function of the positions of the planet and star, we used the technique of 
Winn et al. (2005): we simulated in-transit spectra, and determined the Doppler shifts using 
the same algorithm used on the actual data. The simulations rely on a template spectrum 
(described below) that is meant to mimic the emergent spectrum from a small portion of the 
photosphere. At a given moment of the transit, we denote by e the fractional loss of stellar 
flux, and we denote by Vp the line-of-sight velocity of the occulted portion of the stellar disk. 
To represent the occulted portion of the stellar spectrum, we scaled the template spectrum in 
flux by e and shifted it in velocity by Vp. We subtracted the scaled and shifted spectrum from 
a rotationally-broadened template spectrum and then "measured" the anomalous Doppler 
shift Af. This was repeated for a grid of {e, fp}, and a polynomial function was fitted to 
the resulting grid. We used this polynomial to calculate the anomalous Doppler shift Af 
as a function of e and fp, which are themselves functions of time. Differential rotation was 
ignored, as its effects are expected to be negligible (Gaudi & Winn 2007). 

The template spectrum should be similar to that of TrES-2 but with slightly narrower 
lines because of the lack of rotational broadening. We experimented with two different 
empirical templates based on observations of similar starsj^ finding that both templates gave 
results consistent with the function Af = — e Vp. This function is consistent with the analytic 
expressions of Ohta et al. (2006) and Gimenez (2006), even though those analytic expressions 
do not attempt to account for the spectral deconvolution. It is simpler than the quadratic 
or cubic functions that we have derived for other systems (Winn et al. 2005, 2006, 2007a). 
We do not know the reason for the difference but it is possibly related to the much slower 
projected rotation speed of TrES-2. 

The fitting statistic was 



where /j (obs) and 0"/^- are the flux measurements and uncertainties of Holman et al. (2007), 
and fj(obs) and the radial-velocity measurements and uncertainties from our new 



^The two stars were HD 38858 (Teff = 5726 K, log 5 = 4.51 ± 0.08, [Fe/H] = -0.23 ± 0.04, wsini* = 
0.3 ± 0.5 km s^i) and HD 66428 (Tes = 5752 K, \ogg = 4.49 ± 0.08, [Fe/H] = +0.31 ± 0.04, vsini^ = 
0.0 ± 0.5 km s^^). The stellar parameters are from the SPOCS catalog (Valenti & Fischer 2005). 




(1) 
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data and from O'Donovan et al. (2006). The two model parameters relating to the RM effect 
are the line-of-sight stellar rotation velocity (vsini^), and the angle between the projected 
stellar spin axis and orbit normal (A). The projected spin-orbit angle A ranges from —180° 
to +180°, and is measured counterclockwise on the sky from the projected stellar rotational 
angular- momentum vector to the projected orbital angular-momentum vector (see Ohta et 
al. 2007 or Gaudi & Winn 2007 for a diagram). If we define stellar "north" by the sky 
projection of the stellar angular-momentum vector, then when A = 0° the axes are aligned 
and the planet moves directly "eastward" across the face of the star, for 0° < A < 90° the 
planet moves "northeast," and so forth. 

The other model parameters were the planetary mass (Mp); the stellar and planetary 
radii (i?* and Rp); the orbital inclination (i); the mid-transit time (Tg); and an additive 
constant velocity for each of the two different velocity data sets (71 and 72). We allowed our 
velocities to have a different additive constant from the velocities of O'Donovan et al. (2006) 
in order to account for systematic differences in the spectrograph setup and reduction pro- 
cedures. We fixed the orbital period to be 2.47063 days (Holman et al. 2007). We used a 
Markov Chain Monte Carlo algorithm to solve for the model parameters and their confidence 
limits, with uniform priors on all parameters. This algorithm and our implementation of it 
are described in detail elsewhere (see, e.g., Winn et al. 2007b). The minimum is 1127.6, 
with 1091 degrees of freedom, giving x^/^doi ~ 1.034 and indicating an acceptable fit. 



The RM effect is certainly not obvious in Fig. 1, which shows the entire spectroscopic 
orbit. It is not even very obvious in the middle panel of Fig. 2, which focuses on the 
velocity data around the time of transit. However, our analysis shows that the RM effect 
was indeed detected. As mentioned above, for the best-fitting model, Xmin ~ 1127.6. If the 
parameter vsini^, is set equal to zero, thereby neglecting the RM effect, then x^j^ = 1135.8, 
with the increase of Ax^ = 8.2 arising from the velocity data during the transit. We 
conclude that the RM effect was detected with a signal-to- noise ratio (SNR) of approximately 



V8^ = 2.9. Gaudi & Winn (2007) have given analytic formulas for the signal-to-noise ratio of 



RM observations as a function of the system and telescope parameters, under the assumption 
of Gaussian velocity errors. Using their Eqn. (26) for this case, the forecasted SNR is 2.9, 
in agreement with the actual SNR. 

One might wonder how much this result was influenced by the inclusion of the photo- 
metric data. To check on this, we tried setting aside the photometric data and fltting only 
the 67 radial-velocity data points. We fixed the photometric parameters (Mp, i?*, i, Tc, 
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and P) at the values determined previously. In this case we found xLin = 63.7. If vsini^ 
is set equal to zero, then xLin — 71.9, giving Ax^ = 8.2, just as in the full model fit. This 
confirms that the lowered is an effect of a better fit to the transit velocities, and that the 
uncertainties in the photometric parameters are negligible in this instance. 

The best-fitting model parameters are also consistent with good alignment of the spin 
and the orbit. Specifically, we find A = —9 ± 12 deg, and f sini^ = 1.0 ± 0.6 km s~^, where 
the quoted values are the medians of the a posteriori distributions returned by the MCMC 
algorithm, and the error bars represent 68% confidence limits. Table 2 gives these results, 
along with some other revelant system parameters of TrES-2, for convenience. Visually, the 
RM effect is more apparent in the bottom panel of Fig. 2, in which the orbital velocity 
has been subtracted from the data, and the sampling rate has been effectively doubled by 
inverting the data through the origin {t —t and Af —Av). This works because for 
A ~ 0, the RM waveform is antisymmetric about the origin. 

Figure 3 shows the a posteriori probability distribution for A and the joint distribution 
of A and v sini^. The distribution for A resembles a slightly asymmetric Gaussian function to 
which is added a low-level uniform probability distribution. Although only the region from 
—90° to +90° is shown in Fig. 3, this low- level uniform distribution extends all the way from 
— 180° to +180°. The uniform background corresponds to the very lowest allowed values 
of f sinz^. This makes sense because when the rotation rate is zero, the Rossiter anomaly 
vanishes and A is irrelevant. Values of A between —90° and +90° correspond to prograde 
orbits, for which the stellar and orbital angular momenta are in the same half-plane. The 
integrated probability between —90° and +90° is 98%. We conclude that the TrES-2 orbit is 
prograde with 98% confidence. As an illustration of the constraints provided by our analysis. 
Fig. 4 shows a drawing of the face of the star and the orbit of the transiting planet. 

Our result for v sin i^ is in agreement with the value reported by O'Donovan et al. (2006), 
2.0 ± 1.5 km s~^, which was based on an analysis of the line broadening in an out-of- 
transit spectrum. This finding is also supported by an analysis of our own out-of-transit, 
iodine-free spectra, using the Spectroscopy Made Easy (SME) software package of Valenti & 
Piskunov (1996). The automated analysis gave a formal result of f sinz^ = 0.5 + 0.5 km s~^, 
although the true uncertainty may be larger, since with a disk-integrated spectrum of such 
a slow rotator it is difficult to disentangle the effects of rotation, macroturbulence, microtur- 
bulence, and the instrumental profile. In particular, the SME code assumes "typical" values 
for the turbulent broadening mechanisms that are of the same magnitude as the rotation 



■^We investigated the consequences of accepting the SME resuh at face value, by imposing a Gaussian 
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5. Summciry and Discussion 

We have monitored the apparent Doppler shift of TrES-2 throughout a transit of its giant 
planet and we have detected the Rossiter-McLaughhn effect. Using the available photometric 
and spectroscopic data, we have found good evidence that the orbit is prograde, as are the 
other 6 systems that have been measured (with the possible exception of HD 17156), and as 
are the planets in the Solar system. In this sense, our results for TrES-2 are not surprising. 
However, as mentioned in § 1, some theories of planet migration do predict occasionally large 
misalignments. For example, Nagasawa et al. (2008) investigated a scenario in which a planet 
is scattered into an eccentric, inclined orbit with a small periastron distance (as envisioned 
earlier by Rasio & Ford 1996 and Marzari & Weidenschilling 2002), and subsequently a more 
distant planet forces Kozai oscillations in the inner planet's eccentricity and inclination. If 
the periastron distance is small enough during the high-eccentricity phases, the orbit may 
circularize at a small orbital distance with a substantial inchnation. Nagasawa et al. (2008) 
found that this migration mechanism produces a very broad range of final inclinations, 
including a significant fraction of retrograde orbits. Of course, prograde orbits are also 
permitted in this scenario, and our finding of a prograde orbit for TrES-2 cannot be taken 
as evidence against this mechanism. We raise the issue only to show that a prograde orbit 
was not a foregone conclusion. 

Furthermore, we have shown it is possible to glean this information and measure the 
projected spin-orbit angle to within 12°, even for an 11th magnitude star with a slow pro- 
jected rotation rate. A potentially important application of the RM effect is the detection 
of planets that are too small to be readily detected using other types of ground-based data. 
For example, in many cases of terrestrial planets detected by the Corot or Kepler satellites, 
it will be easier to observe the RM effect than to observe the star's orbital Doppler shift (and 
thereby measure the planet's mass). The theory underlying this idea has been discussed by 
Welsh et al. (2004) and Gaudi & Winn (2007). 

The present work serves to illustrate this point with actual data. If TrES-2 had a 
rotation rate of 5 km s^^ instead of 1 km s^^, but all other stellar and orbital parameters 
were the same, then the quantity and quality of data presented in this paper would permit 
a ~3(7 detection of a planet with a radius ~-\/5 times smaller than TrES-2, or ~6 Earth 
radii. If the transit were equatorial instead of grazing (the best configuration for detecting 
the effect, although not for assessing spin-orbit alignment), the duration of the transit would 



prior constraint on the v sin parameter with mean 0.5 km and standard error 0.5 km s~^. In that case, 
the MCMC analysis gave 68%-confidence ranges of —31 to 1 deg for A and 0.3 to 1.1 km s~^ for vsini^, 
and showed that the orbit is prograde with 95% confidence. The constraint on A was weakened because the 
SME result favors slower rotation rates, for which the sensitivity of the RM waveform to A is reduced. 
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be longer by a factor of ~2 and the amplitude of the RM effect would be larger by a factor 
of ~2, leading to another factor-of-2 improvement in the detectable planet radius (~3 -R©). 
Such a planet would produce a photometric transit depth of only 8 x 10~^, which is smaller 
than the transit depth of any known transiting planet. 

We thank G. Marcy for advice and encouragement, and D. Charbonneau for helpful con- 
versations. We are grateful for support from the NASA Keck PI Data Analysis Fund (JPL 
1326712). We recognize and acknowledge the very significant cultural role and reverence 
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We are most fortunate to have the opportunity to conduct observations from this mountain. 
Access to the Keck telescopes for this project was through the Telescope System Instru- 
mentation Program, and was supported by AURA through the National Science Foundation 
under AURA Cooperative Agreement AST 0132798 as amended. 
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Table 1. Radial Velocities of TrES-2 



HJD 


Radial Velocit [m s~^1 
a la e oci y [m s 


Measurement Uncertainty [m s ^] 


2454216. 


.96699 


42.59 


5. 


.86 


2454216 


.96998 


38.52 


5 


.76 


2454216 


.97930 


37.06 


5 


.43 


2454216 


.98973 


26.69 


5. 


.76 


2454216. 


.99368 


27.35 


5. 


.82 


2454216. 


.99769 


21.46 


5. 


.87 


2454217. 


.00168 


22.74 


6. 


.69 


2454217. 


.00564 


17.82 


5. 


.75 


2454217, 


.00876 


31.44 


6 


.16 


2454217, 


.01102 


16.32 


6 


.06 


2454217, 


.01327 


12.64 


6. 


.24 


2454217, 


.01662 


6.15 


6. 


.05 


2454217, 


.01779 


17.40 


6. 


.13 


2454217, 


.02003 


19.25 


6. 


.35 


2454217, 


.02229 


6.79 


6. 


.24 


2454217, 


.02453 


14.41 


6. 


.35 


2454217, 


.02681 


23.50 


6 


.26 


2454217, 


.02905 


6.40 


6. 


.19 


2454217, 


.03131 


36.58 


6. 


.27 


2454217, 


.03366 


6.90 


6. 


.92 


2454217, 


.03680 


14.89 


6. 


.07 


2454217, 


.03803 


16.89 


6. 


.09 


2454217, 


.04040 


4.04 


6 


.25 


2454217, 


.04266 


7.37 


6 


.17 


2454217, 


.04492 


3.26 


5. 


.91 


2454217, 


.04716 


0.45 


6. 


.44 


2454217, 


.04940 


3.51 


6. 


.27 


2454217, 


.05165 


-0.38 


6. 


.31 


2454217, 


.05403 


2.49 


6. 


.37 


2454217, 


.05648 


7.42 


6 


.13 


2454217, 


.05907 


-8.45 


6 


.16 


2454217, 


.06167 


-4.76 


6. 


.17 


2454217, 


.06426 


-7.03 


6. 


.06 


2454217, 


.06684 


-11.21 


6. 


.06 


2454217, 


.06957 


-8.69 


6 


.32 


2454217, 


.07214 


-6.58 


6 


.15 


2454217, 


.07473 


-17.63 


6 


.14 


2454217, 


.07730 


-19.40 


6 


.14 


2454217, 


.07991 


-25.50 


6. 


.18 


2454217, 


.08260 


-16.69 


6. 


.14 


2454217, 


.08513 


-11.51 


6. 


.33 


2454217, 


.08767 


-16.22 


6 


.21 


2454217, 


.09031 


-28.52 


6 


.52 


2454217, 


.09287 


-18.49 


6 


.62 


2454217, 


.09545 


-9.36 


6 


.50 


2454217, 


.09806 


-25.95 


6. 


.43 


2454217, 


.10066 


-26.17 


6. 


.26 


2454217, 


.10335 


-15.71 


6. 


.36 


2454217, 


.10608 


-15.20 


6 


.34 


2454217, 


.10868 


-21.00 


6 


.26 


2454217, 


.11124 


-30.07 


6 


.31 


2454217, 


.11386 


-22.90 


6 


.24 


2454217, 


.11709 


-35.76 


5. 


.99 


2454217, 


.12110 


-33.87 


6. 


.88 


2454217, 


.12609 


-29.70 


6. 


.79 


2454217, 


.12911 


-26.81 


6. 


.88 



Note. — Column 1 gives the Heliocentric Julian Date at the photon- weighted midex- 
posure time, i.e., weighted by the photon count rate recorded by the HIRES exposure 
meter. 
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Fig. 1. — Radial velocity measurements of TrES-2, from this work and from O'Donovan et 
al. (2006), as a function of orbital phase. The best-fitting values of the systemic velocity 
have been subtracted. The sohd line is the best-fitting model. 

Table 2. System Parameters of TrES-2 



Parameter 


Value ( 


58% Confidence Limits 


References 


P[d] 


2.470621 


±0.000017 


1 


Tc [HJD] 


2453957.63479 


±0.00038 


1 




0.0157 


±0.0003 


1 


b = a cosi/Ri, 


0.8540 


±0.0062 


1 


[Mq] 


0.980 


±0.062 


2 


R* [Rq] 


1.000 


±0.036, -0.033 


2 


Mp [Mjup] 


1.198 


±0.053 


2 


Rp [Rjup] 


1.220 


±0.045, -0.042 


2 


V sin [km s~"^] 


1.0 


±0.6 


This work 


A [deg] 


-9 


±12 


This work 



Note. — (1) Holman et al. (2007); (2) Sozzetti et al. (2007). 
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Fig. 2. — Top. Tlie 2;-band photometry of Holman et al. (2007), averaged into 1.5 min 
bins. The sohd hne is the best-fitting modeL Middle. A close-up of the radial velocity data 
shown in Fig. 1, centered on the midtransit time. Bottom. Same, but the orbital velocity has 
been subtracted and the post-midtransit data {t > 0) have been inverted about the origin 
{t —t and Av — > — A^;), highlighting the Rossiter- McLaughlin anomaly. Filled symbols 
denote data from before midtransit, and open symbols denote data from after midtransit. 
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Fig. 3. — Top. — The probability distribution for A, the angle between the sky projections 
of the orbital axis and the stellar rotation axis. Bottom. — The joint probability distribution 
of A and ?;sini^. The solid dot shows the best-fitting values. The contours represent 68% 
and 95% confidence limits. 
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Fig. 4. — Scale drawing of the TrES-2 system. The relative radii of the bodies and the 
impact parameter of the transit are taken from our best-fitting model. The "north pole" of 
the star is drawn with an arrow, and the curved arc shows the 68%-confidence region for 
its orientation. The angle A is measured clockwise from the projected orbit normal vector 
(vertical dashed line) to the projected stellar north pole (tilted dashed hne). The best-fitting 
value of A is negative. 



